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The finite temperature transition in QCD is studied using Wilson quarks for the cases of Np = 2,3 and 2+1. 
For Np = 2 the transition is smooth in the chiral limit on both TVj = 4 and 6 lattices. For Np = 3, clear two state 
signals are observed for /3 < 4.7 on 8^ X 10 X 4 and 12^ X 4 lattices, which implies the transition is first order for 
rriq < 140 MeV. For 7V_f = 2 + 1 we study two cases of m., ~ 150 and 400 MeV with in.,, = rUri ~ 0. In contrast to 
a previous result with staggered quarks, two state signals are clearly observed for both cases, suggesting a first 
order QCD phase transition in the real world. 



1. INTRODUCTION 

The study of the finite temperature QCD tran- 
sition provides important information for heavy 
ion collisions and the evolution of early Universe: 
Among others the order of the transition is the 
most decisive information. Previous investiga- 
tions show that the order of the chiral transition 
depends sensitively on the number of flavors Np: 
The transition is first order for Np > 3. v^hile 
a continuous transition is suggested for Np = 2. 
Because the s-quark mass is of the same order 
of magnitude as tlie transition temperature, it is 
necessary to include s-quark properly to make a 
prediction for the real world. 

In this report we study the finite temperature 
transition with Wilson quarks. We first present 
the results of the simulations for Np = 2 and 
3 which are the extensions of our previous stud- 
ies[l,2]. We then discuss a more realistic case of 
Np =2 + 1: two degenerate light u and d-quarks 
and a heavier s-quark. Comparison with previous 
results with staggered quarks is also made. We 
mainly perform simulations on 8^ x 10 x N and 
12^ X Nt lattices with Nf (lattice size in the tem- 
poral direction) = 4. Tlie lattice size 10 and 12 in 
one direction is doubled for hadron calculations. 
For Np = 2 we also study the case of Nf = 6. The 
method of simulations is similar to that described 
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2. Np = 2 

For the case of Np = 4 we showed in [2] 
that the finite temperature transition line, Kt- 
line, crosses the cliiral limit line. iTc-line. at 
/3 = PcT — 3.9 - 4.0 and that the transition is 
smooth in the chiral limit. 

In this paper, we extend the study to the case 
of N-t = 6. Our result is shown in Fig. 1. Similar 
to the case of N-t = 4, we find that, when we 
decrease ft toward ftcT — 4.0 - 4.2. ?t?,| decreases 
rapidly and is consistent with zero at /3cT- This 
suggests a continuous transition. 

3. Np = 3 

Appling the same method as in the case of 
Np = 2 to the case of Np =3 at Nt = 4, 
we observed a two-state signal at /3ct — 3.0[2]: 
When we start from a hot configuration, it re- 
mains stable in the deconfining phase. However, 
when we start from a mixed initial configuration, 
the plaquette, the Polyakov loop and decrease 
rapidly with MC-time, showing that the system 
is developing into a confining state. This result 
together with the first order signal observed for 
Np = 6 [1] is consistent with the prediction based 
on universality [3]. 

We here extend the study to finite m^'s by in- 
creasing ft on the iTy-line. Our phase diagram 
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Figure 1. (m^a)^ on the i^c-line for TVt = 4 and 
6 with spatial size 8^ x 10. Tlie locations of the 
crossing point f^cT^^t) are shown by shaded bars. 



thus obtained is given in [4]. We find clear two- 
state signals for /3 = 4.0, 4.5 and 4.7 on both 
8^ X 10 X 4 and 12^ x 4 lattices (Fig. 2). How- 
ever, when we further increase up to /3 > 5.0, we 
find no clear signs of met ast ability. The critical 
value of the quark mass m"'*, up to which we 
get a clear first order signal on these lattices, is 
bounded from below by for /3 = 4.7.^ We find 
TO™*a > 0.175(2) and (m^/TOp)™* > 0.873(6). 

In order to convert the critical quark mass into 
physical units, we use the results aim pa and m^a 
for these /3's. When we plot the masses in the 
confining phase as a function of \jK — l/Kc, 
they are almost independent of /3, Np and Nf for 
/3< 4.7. Fig. 3 shows our results for Np = 2 and 
3. Identifying mp(Kc) with 770 MeV we get a ~ 
0.8 GeV^i for 13 < 4.7. It should be noted that the 
vector meson mass for mq = ~ 150 MeV is 
consistent with = 1020 MeV, indicating that 
our choice of the definition of the quark mass is 



^For niq we use that in the confining phase with the def- 
inition of niq by an axial- vector Ward identity [5,6]. For 
l3 < 5.3. we experience that shows strange behavior in 
the deconhning phase [1]. This can be attributed to an 
effect of 0(a) chiral violation of Wilson fermions. In the 
confining phase, on the other hand, no strange behavior is 
found. A study with an RG improved action shows that 
the masses in the confining phase suffer almost no 0(a) 
effects [7]. 
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Figure 2. Time histories of the plaquette for 
Np = 3 at f3 = 4.7 on a 12^ x 4 lattice. 



close to that of the current quark mass. 

From this value of a we get m"'* > 140 MeV. 

We note that our to"'* is much larger than that 
obtained previously for staggered quarks at Nf = 
4 with spatial 8^ - 16^ lattices [8,9]: two state sig- 



nals are observed for m^a = 0.025, while no clear 
metastabilities are found for TO^a = 0.075. Using 
a result for meson masses at corresponding /3 [10] 
(the values for Np = 4, because they are only 
available) we obtain to"'* ~ 12 - 38 MeV and 
(mT^/mp)"^* ^ 0.42 - 0.58 for staggered quarks. 

4. Np = 2+1 

Let us now turn our attention to the case of 
Np = 2 + 1: rriu = m,d < to^. In order to see what 
happens with realistic quark masses, we study the 
cases vrts — 150 MeV and 400 MeV. Correspond- 
ing values of Kg are obtained from Fig. 3. Keep- 
ing m„ = TOrf ~ we decrease /3 until we hit Kt- 
We find two state signals on an 8^ x 10 x 4 lattice 
for both rris = 150 and 400 MeV and also on a 
12^ X 4 lattice for = 400 MeV (Fig. 4). 

With staggered quarks, Columbia group [9] re- 
ported that no transition occurs at TO„a = m^a = 
0.025, TO^a = 0.1 (to„ = TOd ~ 12 MeV, to^ ~ 50 
MeV using the value of a from [10]). Note in 
this connection that if < m"'* (for degenerate 
Np = 3), we should certainly expect a first order 
transition. Therefore when we recall to"'* > 140 
MeV for Wilson quarks and m"'* 12 - 38 MeV 
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Figure 3. mpa and m^a in the confining phase. 
Open symbols are for Np = 2, /3 = 3.0, 3.5, 4.0, 
4.3, and 4.5 on an 8^ X 10 X 4 lattice. Filled 
symbols are for Np = 3, /3 = 4.0. 4.5 and 4.7 on 
8^ X 10 X 4 and 12^ x 4. Kc{f}) is determined by 
and niq for Np = 2, with errors taking into 
account the difference between two definitions. 



for staggered quarks, the results for Np =2 + 1 
with Wilson and staggered quarks are both con- 
sistent with the results for degenerate Np = 3 
cases, respectively. 

5. CONCLUSION 

We have studied the nature of the finite tem- 
perature transition with Wilson quarks for Np = 
2, 3 and 2+1. For Np = 2 the chiral transition 
is smooth on an TVt = 6 lattice in accord with 
our previous result at N-t — 4. For Np = 3 at 
Nt = 4. clear two state signals are observed on 
the ii'T-hne for rUq < 140 MeV. For iVi? = 2 + 1 
we have studied the cases ~ 150 and 400 MeV 
with m„ = m,i ~ 0. and we have found two state 
signals for both cases. This suggests a first order 
finite temperature transition in the real world. 

The simulations are performed with HITAC 
S820/80 at KEK and with VPP-500 at the Uni- 
versity of Tsukuba. We thank members of KEK 
for their support. We are grateful to C. DeTar 
and K. Rajagopal for useful discussions. This 
work is in part supported by the Grant-in-Aid 
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Figure 4. Time histories of the plaqnette for 
Np = 2 + 1 with TO^ ~ 400 MeV on a 12^ x 4 
lattice. 



of Ministry of Education, Science and Culture 
(N0.O6NPO6OI). 
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